Observation of condensed phases of quasi-planar core-softened colloids 
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We experimentally study the condensed phases of repelling core-softened spheres in two dimen- 
sions. The dipolar pair repulsion between superparamagnetic spheres trapped in a thin cell is 
induced by a transverse magnetic field and softened by suitably adjusting the cell thickness. We 
scan a broad density range and we materialize a large part of the theoretically predicted phases in 
systems of core-softened particles, including expanded and close-packed hexagonal, square, chain- 
like, stripe/labyrinthine, and honeycomb phase. Further insight into their structure is provided by 
Monte Carlo simulations. 
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One of the ageless challenges of material science is 
to engineer structure formation. Recently, of particu- 
lar interest are the various crystalline phases of colloidal 
particles [T] which may be used for self-assembled nano- 
materials |2, photonic crystals [5], macroporous poly- 
mer membranes jj], biochemical applications [5], etc. In 
this context, it is imperative to understand the link be- 
tween the interparticle potential and the phase behavior 
of the ensemble. Experimentally, this link can be stud- 
ied readily using table-top equipment which provides the 
full microscopic real-time structural information, where- 
from the subtle relations between the microscopic and 
the macroscopic properties of the system can be recon- 
structed. Complementary insight can be obtained by in- 
verse methodology to identify the optimal pair potential 
that produces a given target structure [6J. 

A basic limitation of these efforts is the spherical shape 
typical for many colloids and the ensuing isotropy of the 
interparticle potential. This is the main reason why the 
most common ordered structures are close-packed, i.e., 
the hexagonal and the face-centered cubic lattice [71 E] 
in two and three dimensions, respectively. Conceivably, 
a richer phase diagram can be induced by anisotropic in- 
teractions due to, e.g., surface treatment of particles [9J, 
external fields [10] or a liquid-crystalline solvent [11]. An- 
other route to the more open lattices are isotropic pair 
interactions with a radial profile characterized by two 
length scales, such as the combination of hard-core and 
penetrable-sphere repulsion [12]. If the shoulder/core di- 
ameter ratio exceeds about 2, this pair potential sta- 
bilizes a range of mesophases intervening between the 
fluid and the close-packed crystal. In two dimensions, 
the theoretical phase sequence includes loose- and close- 
packed hexagonal lattice; monomer, dimer, and trimer 
fluids; stripe and labyrinthine phases; honeycomb lat- 
tice, etc. [12 HSl [H]- Very similar behavior has been 
predicted numerically in paramagnetic particles confined 
to a plane and interacting with a dipolar repulsion in- 
duced by a transverse magnetic field and softened by a 
Lcnnard- Jones interaction [TSj. For large shoulder/core 
diameter ratios, the set of mesophases reduces to micel- 
lar, lamellar, and inverted micellar structure [TB] . 



The theoretical insight into the mesophases formed by 
particles with a core-softened isotropic repulsive pair po- 
tential in two dimensions is reasonably comprehensive 
but the physical evidence of their stability is still lack- 
ing. In this Letter, we study the phase sequence of such a 
system experimentally using superparamagnetic spheres, 
tailoring the profile of interaction with external magnetic 
field and spatial constraints. We discover a host of the 
predicted structures, and we corroborate the observations 
with numerical simulations. 

We used a thin wedge-shaped cell formed by nearly 
parallel glass plates and filled by a water solution of 
1.05 /im superparamagnetic spheres (Dynabeads, My- 
One Carboxy). A small (5 mg/ml) amount of SDS was 
added to suppress the van der Waals interactions and 
prevent sticking of the colloids. Like in related experi- 
ments [T71 [IHl Uni Uni 121] , the interparticle repulsion was 
induced by external magnetic field. The experimental 
setup is comprised of laser tweezers system built around 
an inverted optical microscope (Zeiss Axiovert 200M) and 
a CW diode pumped Nd:YAG laser (Coherent, Com- 
pass 2500MN) , and equipped with three orthogonal pairs 
of Helmholtz coils driven by a computer-controlled cur- 
rent source. Multiple laser traps used to manipulate 
the spheres are created by time-sharing using acousto- 
optic defiectors (IntraAction, DTD-274HA6) driven by 
beam-steering controller (Aresis, Tweez). The cell is il- 
luminated with a halogen lamp and the micrographs are 
recorded in the bright field with a fast CMOS camera 
(PixeLINK, PL-A741). Image analysis is performed of- 
fline with custom-made particle tracking software. 

The system's key feature which softens the repulsion 
between the induced magnetic dipoles of the spheres is 
the cell thickness: If it is somewhat larger than the sphere 
diameter, the centers of two near-by spheres are not re- 
stricted to a plane perpendicular to the magnetic field. 
This makes their interaction at small separations less re- 
pulsive (or even attractive) compared to spheres lying in 
the plane. The pair potential is 
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where K = ttct^x^B'^ /XAAfiQ is the interaction constant 
which depends on the magnetic induction B and the mag- 
netic susceptibihty x of the particles; r and z are the in- 
plane and the vertical separations of the spheres, respec- 
tively, and cr is their diameter (Fig.jl]). It is convenient to 
write the cell thickness as a + h: Thus h measures the de- 
viation from a truly planar system. In the thin part of the 
cell where h ^ 0, U reduces to K/r^, whereas for h > 
the repulsion is reduced at small distances by the relative 
vertical shift of the spheres. The in-plane interparticle 
force is Fr = -dU/dr = ~3Kr {Az^ - r^) / {r^ + z^Y^^ . 
For spheres in contact in a cell of thickness a + h < 2(t, 
Fr ~ —3K\/a'^ — h"^ (bh? — u^) /cr^ which is attractive 

for h > hm = cr/\/5 ~ 0.447cr. The softening of the re- 
pulsion is most pronounced at cell thicknesses somewhat 
smaller than a + hm ■ 

To verify that the in-plane pair interaction between 
the spheres is indeed described by the above law, we 
measured it using two isolated spheres, one attached to 
the lower glass plate and the other one (probe) confined 
laterally by a weak optical trap and pushed against the 
upper plate by the light pressure so that their vertical 
separation was fixed to h. The probe was repeatedly 
slowly dragged back and forth from the fixed sphere, and 
its trajectories were recorded and averaged to accurately 
determine the displacement of the probe particle from 
the trap center as a function of interparticle distance. 
To calculate the magnetic force on the probe from the 
displacement profile, the stiffness of the laser potential 
was determined by analyzing the Brownian motion of 
the probe sphere with the magnetic field turned off. The 
procedure was repeated at a range of cell thicknesses and 
the three representative force profiles shown in Fig. [l] 
demonstrate that the system behaves as expected. 

Having confirmed that the pair potential is given by 
Eq. ([!]), we studied the phase behavior of the system 
across a broad density range. We filled the cell with a 
dense suspension of spheres and we first located the mea- 
suring site with cell thickness of about a + hm- As the 
cell thickness cannot be measured directly, we used the 
following procedure. In absence of the magnetic field, we 
herded the spheres at a given location into contact such 
that they formed clusters. We then turned on the mag- 
netic field which caused the clusters to disintegrate into 
smaller chunks. In the thick part of the cell where the 
potential was attractive at small separations the cluster 
disintegration was partial, whereas in the thin part with 
purely repulsive interactions the clusters disintegrated 
completely. By scanning the region we identified the site 
with the desired cell thickness. 

At this site, we increased the sphere density n in the 
part of the cell covered by camera's field of view using 
the laser tweezers operating at high power to locally heat 
the suspension. The heating induced a hydrodynamic 
flow which dragged a large number of spheres (typically 
over 10^) towards the trap, thereby increasing the local 
sphere density almost to close-packing. Then the laser 
was turned off which stopped the flow, and the magnetic 




FIG. 1: Measured variability of the in-plane pair force profiles: 
r^'* repulsion (thin cell, h <C hm', circles), softened repulsion 
{h ~ hm', squares), and oversoftened interaction (thick cell, 
h > hm', triangles) with the attractive part at small separa- 
tions. For h = hm, the in-plane hard-core diameter of the 
spheres is ~ 0.89 /^m. The small-separation part of the force 
profile is not accessible with our method as the laser would 
act on both spheres, rendering the measurement impossible. 
In our experiment, B is typically around 10 mT and the in- 
duced pair forces are of the order of 0.1 pN at separations 
comparable to sphere diameter. — Solid lines are theoretical 
fits for h = 0, h = 0.46 /im and h — 0.72 fim, respectively. 
Inset: Experimental geometry in cross-section. 



field was turned on. The high-density region underwent 
an adiabatic expansion so the system could be studied at 
decreasingly lower colloidal densities. A sizable change 
of the density took place on a time scale of about 10 s. 
This is much longer than the typical diffusion time of 
particles < 10 ms [22] so that the expansion was slow 
enough to ensure quasi-equilibrium at all times, and the 
hydrodynamic interactions due to expansion ~ lO^'^fc^T 
are negligible compared to the magnetic repulsion. 

Micrographs of the most interesting phases observed 
at various 2D volume fractions i] = iva^n/A shown in 
Fig. |2] include i) fluid phase {-q — 0.01); ii) expanded 
hexagonal lattice (77 — 0.12); iii) coexisting expanded 
hexagonal and square lattice (77 — 0.23); iv) coexist- 
ing expanded square lattice and chain phase (77 = 0.31); 
v) chain phase with locally aligned flnite-length strings 
(ry = 0.34); vi) stripe/labyrinthine structure formed pre- 
dominantly by a single cluster of interlaced strings of 
touching spheres ij] — 0.39); and vii) coexisting hon- 
eycomb and dense square lattice at very high volume 
fraction (77 — 0.54). The observed mesophases are re- 
markably close to those found in the numerical simula- 
tions reported in Ref. IF] although the pair interaction 
is not exactly the same; a much more idealized hard- 
core/soft-shoulder interaction also gives a similar phase 
sequence [12l [14] . This suggests that the mechanisms at 
work as well as the structures they produce are rather 
robust. 



3 




FIG. 2: Micrographs of the representative mesophases induced by varying the density (top row; B = 12.5 mT which gives 
k ~ 400) and the simulation snapshots (bottom row) labelled by the corresponding 2D volume fractions of colloids. Due to 
optical image distortion, the colloids on the micrographs appear slightly larger than their true size; in some of the micrographs, 
patches of the underlying lattices are emphasized to guide the eye. — The reduced interaction strength of fc = 375 was used in 
simulations at all but the largest two volume fractions where k was set to 125 to improve convergence; the number of particles 
iV = 1000 except at 77 = 0.60 where TV = 200. 



Nonetheless, our system is not perfectly two- 
dimensional, and the vertical position of a sphere is a 
restricted variable which does affect the symmetry of the 
in-plane equilibrium configuration. To understand the 
phase sequence in more detail, we performed 3D Monte 
Carlo simulations of up to TV = 1000 spheres with the 
dipole-dipole pair repulsion in a cell of thickness a + hm- 
We varied the volume fraction and we focused on the 
low temperature limit corresponding to large reduced in- 
teraction strengths k — K/ksTa^; in the experiment, 
k « 400. We used periodic boundary conditions in sev- 
eral simulation box geometries; after reaching equilib- 
rium (typically in a few million steps) 50 x 10^ averaging 
steps were performed to evaluate the energy per particle, 
the radial distribution function, and the static structure 
factor. The selected simulation snapshots in Fig. [2] are 
chosen to best reproduce the phase sequence seen in the 
micrographs. The agreement is very good, which also 
demonstrates that any effect of polydispersity of the su- 
perparamagnetic spheres used is unimportant. 

The insight provided by the simulations goes beyond 
the experimental micrographs: It includes the vertical 
positions of the spheres, which are uncorrelated and 
evenly distributed across the available range as long as 
the volume fraction does not exceed about 0.05. How- 
ever, for rj > 0.1 the distribution becomes bimodal 
with pronounced peaks corresponding to spheres touch- 
ing the bottom and the top plate, and the system re- 
sembles an off-lattice two-state spin ensemble with dom- 
inant nearest-neighbor antiferromagnetic interactions. 
The nearest-neighbor interactions are frustrated by the 
hexagonal lattice but the stripe, honeycomb, and square 
lattice with 2, 3, and 4 regularly arranged nearest neigh- 
bors, respectively, are compatible with alternating up- 
down positions of spheres. We calculated the energies of 
the stripe, square, and hexagonal phase by fixing their 



in- plane structure [23J. The ground states of the non- 
frustrated phases, such as the checkerboard structure of 
the square lattice, are ordered and unique, and their en- 
ergies were calculated using a lattice sum, whereas the 
energy of the hexagonal lattice was evaluated by anneal- 
ing the vertical positions numerically in systems of up to 
TV = 5000 spheres. 

The results are shown in Fig.|3] At small 77, the hexag- 
onal phase has the lowest energy per particle among the 
three candidate phases, the square phase is stable at in- 
termediate volume fractions, and at largest volume frac- 
tions considered the stripes are energetically most favor- 
able. In Fig. [3) we also plotted the energies of the equi- 
librium states obtained by the full Monte Carlo analysis 
that produced the snapshots shown in Fig. [2j Up to 
T] 0.4, the agreement of the energies calculated us- 
ing lattice sums and simulations, respectively, is quite 
good. In this regime, the simulations as well as the ex- 
periments usually produce coexistence of the hexagonal 
and the square lattice rather than a pure lattice (Fig. |2| . 
This can be understood based on the lattice sum results 
which predict a very similar dependence of the energies 
of the two lattices on 77, and thus the coexistence regime 
obtained by Maxwell construction should extend over 
a broad density range. — At volume fractions beyond 
w 0.4, the difference between the energies of the model 
structures and those obtained by simulations slightly ex- 
ceeds the error bar of the latter suggesting that the stripe 
phase with perfectly parallel straight stripes is probably 
not the ground state and many turns and junctions may 
be an equilibrium feature of the system. 

In conclusion, we report an experimental study where 
a core-softened repulsive interaction was induced in 
micrometer-size colloidal particles using the external 
magnetic field and fine-tuned by controlling the spatial 
constraints. Depending on the density, the system forms 
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several self-assembled mesophases — the square, hexag- 
onal, and honeycomb lattices as well as the labyrinthine 
structure — , thereby experimentally validating the the- 
oretical predictions pertaining to similar pair interac- 
tions [m [131 Hn] . The observed two-dimensional struc- 
tures could be used for the fabrication of templates to 
promote the growth of colloidal crystals [24} 125]. 
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FIG. 3: Energy per particle of the ground states of hexago- 
nal (broad shaded line; the line thickness represents the nu- 
merical inaccuracy estimated by varying both the number of 
particles and the coupling constant), square (solid line), and 
stripe phase (dashed line) compared to the energies obtained 
by the full Monte Carlo simulation (open circles with error 
bars). Also shown are the snapshots for = 0.10 (hexagonal 
lattice), ri — 0.30 (square lattice), and rj = 0.52 (labyrinthine 
structure). The spheres' vertical positions are encoded by 
shades of gray (light gray: z = 0, black: z — h). The checker- 
board pattern of the square lattice is clearly visible, and so is 
the alternating up-down intra-stripe order of the labyrinthine 
structure. In the hexagonal lattice, the vertical positions of 
the spheres are disordered. 



We thank H. H. von Griinberg for stimulating dis- 
cussions and B. Kavcic for technical assistance. This 
work was supported by the European Commission (J. 
D., Grant MERG-031089 and by the Slovenian Research 
Agency (D. B. and I. P., Grant Jl-6502; J. D. and P. Z., 
Grant Pl-0055). 



[1] Y. Xia, B. Gates, Y. Yin, and Y. Lu, Adv. Mater. 12, 
693 (2000). 

[2] J. H. Fendler, Chem. Mater, 8, 1616 (1996). 

[3] K. Busch and S. John, Phys. Rev. E 58, 3896 (1998). 

[4] P. Jiang, K. S. Hwang, D. M. Mittleman, J. F. Bertone, 

and V. L. Colvin, J. Am. Chem. Soc. 121, 11630 (1999). 
[5] H. Kawaguchi, Prog. Polym. Sci. 25, 1171 (2000). 
[6] M. C. Rechstman, F. H. Stillinger, and S. Torquato, 

Phys. Rev. Lett. 95, 228301 (2005). 
[7] K. Kremer, M. O. Robbins, and G. S. Grest, Phys. Rev. 

Lett. 57, 2694 (1986). 
[8] S. C. Mau and D. A. Huse, Phys. Rev. E 59, 4396 (1999). 
[9] D. R. Nelson, Nano Lett. 2, 1125 (2002). 
[10] A. Yethiraj and A. van Blaaderen, Nature 421, 513 

(2003). 

[11] H. Stark, Phys. Rep. 351, 387 (2001). 
[12] G. Malescio and G. Pellicane, Nat. Mater. 2, 97 (2003). 
[13] Y. Norizoe and T. Kawakatsu, Europhys. Lett. 72, 583 
(2005). 

[14] E. A. Jagla, J. Chem. Phys. 110, 451 (1999). 
[15] P. J. Camp, Phys. Rev. E 68, 061506 (2003). 
[16] M. Glaser, G. M. Grason, R. D. Kamien, A. Kosmrlj, 

C. D. Santangelo, and P. Ziherl, Europhys. Lett. 78, 

46004 (2007). 

[17] K. Zahn, R. Lenke, and G. Maret, Phys. Rev. Lett. 82, 
2721 (1999). 

[18] R. Bubeck, S. Neser, C. Bechinger, and P. Leiderer, Prog. 

Colloid Polym. Sci. 110, 41 (1998). 
[19] W. Wen, L. Zhang, and P. Sheng, Phys. Rev. Lett. 85, 



5464 (2000). 

[20] C. Eisenmann, U. Gasser, P. Keim, and G. Maret, Phys. 
Rev. Lett. 93, 105702 (2004). 

[21] A. Erbe and L. Baraban, private communication. 

[22] In equilibrium, a thermally excited particle in the po- 
tential well formed by its neighbors fluctuates around 
its mean position by about d ~ [ksTa^ f{ri)/38AKY^'^ 
where /(t;) is a numerical factor which depends on the 
volume fraction and equals 1 at close-packing and ~ 1000 
at 77 « 0.1. The typical time needed to traverse this dis- 
tance is T ~ /D, where D = kBT/Sivua is the diffusion 
constant; 1/ is the viscosity. In our case, r ~ 10 /is at 
77 = 0.5 and r ~ 10 ms at = 0.1. 

[23] We analyzed all stripe structures of width 1 by varying 
the intra-stripe separation, the stripe spacing, and the 
staggering of neighboring stripes. At any volume fraction, 
the minimal-energy configuration consists of close-packed 
stripes, and the neighboring stripes are out of registry by 
1 sphere diameter such that the stripe structure resem- 
bles a square lattice compressed along one lattice vector 
and stretched along the other one as much as possible. 

[24] A. van Blaaderen, R. Ruel, and P. Wiltzius, Nature 385, 
321 (1997). 

[25] K. H. Lin, J. C. Crocker, V. Prasad, A. Schofield, 
D. A. Weitz, T. C. Lubensky and A. G. Yodh, Phys. 
Rev. Lett. 85, 1770 (2000). 



